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Abstract 

This paper discusses hadron energy reconstruction for the ATLAS 
barrel prototype combined calorimeter (consisting of a lead-liquid ar- 
gon electromagnetic part and an iron-scintillator hadronic part) in the 
framework of the non-parametrical method. The non-parametrical me- 
thod utilizes only the known e/h ratios and the electron calibration con- 
stants and does not require the determination of any parameters by a min- 
imization technique. Thus, this technique lends itself to an easy use in a 
first level trigger. The reconstructed mean values of the hadron energies 
are within ±1% of the true values and the fractional energy resolution is 
[(58 ± 3)%/VE + (2.5 ± 0.3)%] © (1.7 ± 0.2) /E. The value of the 
e/h ratio obtained for the electromagnetic compartment of the combined 
calorimeter is 1.74 ± 0.04 and agrees with the prediction that e/h > 1.7 
for this electromagnetic calorimeter. Results of a study of the longitudi- 
nal hadronic shower development are also presented. The data have been 
taken in the H8 beam line of the CERN SPS using pions of energies from 
10 to 300 GeV. 

Codes PACS: 29.40. Vj, 29.40.Mc, 29.85,+c. 

Keywords: Calorimetry, Combined Calorimeter, Shower Counter, Com- 
pensation, Energy Measurement, Computer Data Analysis. 
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1 Introduction 



The key question for calorimetry in general, and hadronic calorimetry in par- 
ticular, is that of energy reconstruction. This question becomes especially 
important when a hadronic calorimeter has a complex structure incorporating 
electromagnetic and hadronic compartments with different technologies. This 
is the case for the central (barrel) calorimetry of the ATLAS detector which 
has the electromagnetic liquid argon accordion and hadronic iron- scintillator 
Tile calorimeters [[]], 0, |3|] . A view of the ATLAS detector, including the two 
calorimeters, is shown in Fig. []]. 

In this paper, we describe a non-parametrical method of energy reconstruc- 
tion for a combined calorimeter known as the e/h method, and demonstrate 
its performance using the test beam data from the ATLAS combined prototype 
calorimeter. For the energy reconstruction and description of the longitudinal 
development of a hadronic shower, it is necessary to know the e/h ratios, the 
degree of non-compensation, of these calorimeters. Detailed information about 
the e/h ratio for the ATLAS Tile barrel calorimeter is presented in [Q, |], |5], [7p 
while much less was done so far for the liquid argon electromagnetic calorime- 
ter |9], [TTJ. An additional aim of the present work, then, is to also determine 
the value of the e/h ratio for the electromagnetic compartment. 

Another important question for hadron calorimetry is that relating to the 
longitudinal development of hadronic showers. This question is especially im- 
portant for a combined calorimeter because of the different degrees of non- 
compensation for the separate calorimeter compartments. Information about 
the longitudinal hadronic shower development is very important for fast and 
full hadronic shower simulations and for fast energy reconstruction in a first 
level trigger. This work is also devoted to the study of the longitudinal hadronic 
shower development in the ATLAS combined calorimeter. 

This work has been performed using the 1996 combined test beam data 
JTT] , [12| ] taken in the H8 beam line of the CERN SPS using pions of energies 
from 10 to 300 GeV. 

2 Combined Calorimeter 

The combined calorimeter prototype setup is shown in Fig. |2|, along with a def- 
inition of the coordinate system used for the test beam. The LAr calorimeter 
prototype is housed inside a cryostat with the hadronic Tile calorimeter proto- 
type located downstream. 
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The beam line is in the YZ plane at 12 degrees from the Z axis. With this 
angle the two calorimeters have an active thickness of 10.3 interaction lengths 
(A/). The beam quality and geometry were monitored with a set of scintilla- 
tion counters SI - S4, beam wire chambers BC1 - BC3 and trigger hodoscopes 
(midsampler) placed downstream of the cryostat. To detect punchthrough parti- 
cles and to measure the effect of longitudinal leakage a "muon wall" consisting 
of 10 scintillator counters (each 2 cm thick) was located behind the calorimeters 
at a distance of about 1 metre. 

The liquid argon electromagnetic calorimeter prototype consists of a stack 
of three azimuthal modules, each module spanning 9° in azimuth and extending 
over 2000 mm along the Y direction. The calorimeter structure is defined by 
2.2 mm thick steel-plated lead absorbers folded into an accordion shape and 
separated by 3.8 mm gaps filled with liquid argon. The signals are collected 
by three-layer copper-polyamide electrodes located in the gaps. The calorime- 
ter extends from an inner radius of 1315 mm to an outer radius of 1826 mm, 
representing (in the Z direction) a total of 25 radiation lengths (X ), or 1.22 
Xi for protons. The calorimeter is longitudinally segmented into three com- 
partments of 9 X , 9 X and 7 X , respectively. The 77 x <\> segmentation is 
0.018 x 0.02 for the first two longitudinal compartments and 0.036 x 0.02 for 
the last compartment. Each read-out cell has full projective geometry in r\ and 
<p. The cryostat has a cylindrical shape, with a 2000 mm internal diameter (filled 
with liquid argon), and consists of an 8 mm thick inner stainless- steel vessel, 
isolated by 300 mm of low-density foam (Rohacell), which is itself covered by 
a 1.2 mm thick aluminum outer wall. A presampler was mounted in front of 
the electromagnetic calorimeter. The presampler has fine strips in the rj direc- 
tion and covers »llx8inr]X0 LAr calorimeter cells in the region of the 
beam impact. The active depth of liquid argon in the presampler was 10 mm 
and the strip spacing 3.9 mm. Early showers in the liquid argon were kept to 
a minimum by placing light foam material (Rohacell) in the cryostat upstream 
of the LAr electromagnetic calorimeter. The total amount of material between 
BC3 and LAr calorimeter is near 0.2Aj. More details about this prototype can 
be found in [0, 

The hadronic Tile calorimeter is a sampling device which uses steel as the 
absorber and scintillating tiles as the active material [@]. A conceptual design 
of this calorimeter geometry is shown in Fig. ^. The innovative feature of the 
design is the orientation of the tiles which are placed in planes perpendicular 
to the Y direction JT3[]. The absorber structure is a laminate of steel plates 
of various dimensions stacked along Y. The basic geometrical element of the 
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stack is denoted as a period. A period consists of a set of two master plates 
(large trapezoidal steel plates, 5 mm thick, spanning along the entire Z dimen- 
sion) and one set of spacer plates (small trapezoidal steel plates, 4 mm thick, 
100 mm wide along Z). During construction, the half-period elements are pre- 
assembled starting from an individual master plate and the corresponding 9 
spacer plates. The relative position of the spacer plates in the two half periods 
is staggered in the Z direction, to provide pockets in the structure for the sub- 
sequent insertion of the scintillating tiles. Each stack, termed a module, spans 
27r/64 in the azimuthal angle (X dimension), 1000 mm in the Y direction and 
1800 mm in the Z direction (about 9 A/ or about 80 X ). The module front 
face, exposed to the beam particles, covers 1000x200 mm 2 . The scintillating 
tiles are made out of polystyrene material of thickness 3 mm, doped with scin- 
tillating and wavelength- shifting dyes. The iron to scintillator ratio is 4.67 : 1 
by volume. The tile calorimeter thickness along the beam direction at the in- 
cidence angle of 12° (the angle between the incident particle direction and the 
normal to the calorimeter front face) corresponds to 1.5 m of iron equivalent 
length. 

Wavelength shifting fibers collect the scintillation light from the tiles at both 
of their open (azimuthal) edges and transport it to photo-multipliers (PMTs) at 
the periphery of the calorimeter (Fig. Q). Each PMT views a specific group of 
tiles through the corresponding bundle of fibers. The prototype Tile calorimeter 
used for this study is composed of five modules stacked in the X direction, as 
shown in Fig. |2|. 

The modules are longitudinally segmented (along Z) into four depth seg- 
ments. The readout cells have a lateral dimension of 200 mm along Y, and 
longitudinal dimensions of 300, 400, 500, 600 mm for depth segments 1-4, 
corresponding to 1.5, 2, 2.5 and 3 \j, respectively. Along the X direction, the 
cell sizes vary between about 200 and 370 mm depending on the Z coordinate 
(Fig. @). More details of this prototype can be found in [Jl], [14], |3| g, 0, [Hp. 
The energy release in 100 different cells was recorded for each event [|I4|]. 

The data have been taken in the H8 beam line of the CERN SPS using 
pions of energy 10, 20, 40, 50, 80, 100, 150 and 300 GeV. We have applied 
some cuts similar to [O, [T^p in order to eliminate the non-single track pion 
events, the beam halo, the events with an interaction before the liquid argon 
calorimeter, and the electron and muon events. The set of cuts adopted is as 
follows: single-track pion events were selected by requiring the pulse height 
of the beam scintillation counters and the energy released in the presampler of 
the electromagnetic calorimeter to be compatible with that for a single particle; 
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the beam halo events were removed with appropriate cuts on the horizontal and 
vertical positions of the incoming track impact point and the space angle with 
respect to the beam axis as measured with the beam chambers; a cut on the total 
energy rejects incoming muons. 



3 The e/h Method of Energy Reconstruction 

An hadronic shower in a calorimeter can be seen as an overlap of a pure elec- 
tromagnetic and a pure hadronic component. In this case an incident hadron 
energy is E = E e + Eh- The calorimeter response, R, to these two components 



is usually different [[18], [190 and can be written as: 

R = e-E e + h-E h , (1) 

where e (h) is a coefficient to rescale the electromagnetic (hadronic) energy 
content to the calorimeter response. A fraction of an electromagnetic energy of 
a hadronic shower is f n o = E e /E, than R = e ■ f w o ■ E + h ■ (E — f w o ■ E) = 
e- [1 + (e/h—l)-f 7T o]/(e/h)-E. From this one can gets formulae for an incident 
energy 

*-;■&)•*• (2) 

where 

e/h 



— (3) 
nj 1 + (e/h - 1) • fa 

The dependence of f n Q from the incident hadron energy can be parameterized 
as inRef. [£Dp: 

f 7T o = k-\nE. (4) 

In the case of the combined setup described in this paper, the total energy 
is reconstructed as the sum of the energy deposit in the electromagnetic com- 
partment (ElAtX me deposit in the hadronic calorimeter (Etiu), and that in the 
passive material between the LAr and Tile calorimeters (E dm ). Expression @ 
can then be rewritten as: 

E = E L A r + E dm + E T u e = ( — ) Rlat + E dm H ( — 1 Rnu 



eLAr\Kj LAr . e T ile \TT / Tile 

(5) 

where Rlat (Rrue) is the measured response of the LAr (Tile) calorimeter com- 
partment and 1 / exile an d 1 / tLAr are energy calibration constants for the LAr 
and Tile calorimeters respectively [ |TT| ] . 
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Similarly to the procedure in Refs. [(TTJ, ^Tp, the Ed m term, which accounts 
for the energy loss in the dead material between the LAr and Tile calorimeters, 
is taken to be proportional to the geometrical mean of the energy released in 
the third depth of the electromagnetic compartment and the first depth of the 
hadronic compartment {Edm = ol ■ ^Elav,3 • -EriZe.i)- The validity of this ap- 
proximation has been tested using a Monte Carlo simulation along with a study 
of the correlation between the energy released in the midsampler and the E dm 

ng,E2rai. 



The ratio (e/h)Tu e = 1-30 ± 0.03 has been measured in a stand-alone test 
beam run [g| and is used to determine the (e/ir)Tu e term in equation |5]. To 
determine the value of the 1/erue constant we selected events which started 
showering only in the hadronic compartment, requiring that the energy de- 
posited in each sampling of the LAr calorimeter and in the midsampler is 
compatible with that of a single minimum ionization particle. The result is 
l/e Til e = 0.145 ±0.002. 

The response of the LAr calorimeter has already been calibrated to the elec- 
tromagnetic scale; thus the constant 1/e LAr = 1 [ |TT| , |T2"| ] . The value of(e/ti) LAr 
has been evaluated using the data from this beam test, selecting events with well 
developed hadronic showers in the electromagnetic calorimeter, i.e. events with 
more than 10% of the beam energy in the electromagnetic calorimeter. Using 
the expression (Q), the (e/n) LAr ratio can be written as: 

e \ E beam — E dm — E Ti i e 

= 5 7 • W 

11 / LAr t^LArl ^LAr 

Fig. [| shows the distributions of the [ejix^LAr ratio for different energies, and 
the mean values of these distributions are plotted in Fig. |] as a function of the 
beam energy. From a fit to this distribution using expression (|J) and (Q) we 
obtain (e/h) LAr = 1.74 ±0.04 and k = 0.108 ±0.004, thereby taking (e/h) LAr 
to be energy independent. For a fixed value of the parameter k = 0.11 [j20|], 
the result is (e/h) LAr = 1-77 ± 0.02. The quoted errors are the statistical ones 
obtained from the fit. The systematic error on the (e/h) LAr ratio, which is a 
consequence of the uncertainties in the input constants used in the equation 
@ as well as of the shower development selection criteria, is estimated to be 
±0.04. 

Figure |^ compares our values of the (e /if) LAr ratio to the ones obtained 



in Refs. [5| |T0p using a weighting method. The results are in good agree- 
ment below 100 GeV but disagree above this energy because the weighting 
method leads to a distortion of the (e/7r) Lj4r ratios. Despite this disagreement, 
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fitting expression @ to the old data leads to (e/h) em = 1.73 ± 0.10 for [Q] and 
(e/ h) em = 1.64 ± 0.18 for [|10p (parameter k fixed at 0.1 1). These values are in 
agreement with our result within error bars. 

In the Ref. fl2"0] ] it was demonstrated that the e/h ratio for non-uranium 
calorimeters with high-Z absorber material is satisfactorily described by the 
formula: 

e _ e/mip 
h 0.41 + f n ■ n/mip ' 
where f n is a constant determined by the Z of the absorber (for lead /„ = 0.12) 



[ |24| , [250 , and e/mip and n/mip represent the calorimeter response to electro- 
magnetic showers and to MeV-type neutrons, respectively. These responses are 
normalized to the one for minimum ionizing particles. The Monte Carlo cal- 
culated e/mip and n/mip values [ |T8"| ] for the lead liquid argon electromagnetic 
calorimeter [ |2~6| ] are e/mip = 0.78 and n/mip < 0.5, leading to e/h > 1.66. 
The measured value of the (e/h) em ratio agrees with this prediction. Using 
expression (J7|) and measured value of e/h we can find that n/mip is ~ 0.3. 

Formula (|7|) indicates that e/mip is very important for understanding com- 
pensation in lead liquid argon calorimeters. The degree of non-compensation 
increases when the sampling frequency is also increased [ [241 ] ■ A large fraction 
of the electromagnetic energy is deposited through very soft electrons (E < 1 
MeV) produced by Compton scattering or the photoelectric effect. The cross 
sections for these processes strongly depend on Z and practically all these pho- 
ton conversions occur in the absorber material. The range of the electrons pro- 
duced in these processes is very short, ~ 0.7 mm for 1 MeV electron in lead. 
Such electrons only contribute to the calorimeter signal if they are produced 
near the boundary between the lead and the active material. If the absorber 
material is made thinner this effective boundary layer becomes a larger fraction 
of the total absorber mass and the calorimeter response goes up. This effect 
was predicted by EGS3 simulation p7[]. It leads to predictions for the GEM 
p8| ] accordion electromagnetic calorimeter (1 mm lead and 2 mm liquid argon) 
that e/mip = 0.86 and e/h > 1.83. The Monte Carlo calculations also pre- 
dict that the electromagnetic response for liquid argon calorimeters (due to the 
larger Z value of argon) is consistently larger than for calorimeters with plastic - 
scintillator readout. The signal from neutrons (n/mip) is suppressed by a factor 
0.12 and the n — p elastic scattering products do not contribute to the signal of 
liquid argon calorimeters. These detectors only observe the 7's produced by 
inelastic neutron scattering (thermal neutron capture escapes detection because 
of fast signal shaping) [E3p . 
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To use expression (|5|) for reconstructing incident hadron energies, it is nec- 
essary to know the (e/ir)Tiie an d (e/^LAr ratios, which themselves depend on 
the hadron energy. For this purpose, a two cycle iteration procedure has been 
developed. In the first cycle, the {e/i^TUe rat i° is iteratively evaluated using the 
expression: 

( e \ = (e/h) Ti ie (g) 

V 71 "/ Tile 1 + (i e / h )Tile - 1) • A; • In (l/erife • {e/^Tile ■ Rrue) ' 

using the value of (e/ir) Ti i e from a previous iteration. To start this procedure, a 
value of 1.13 (corresponding to /° = 0.11 ln(100 GeV)) has been used. 

In the second cycle, the first approximation of the energy, E, is calculated 
using the equation (|5|) with the (e/ir)Tu e ratio obtained in the first cycle and 
the [ej^LAr ratio from equation (Bl), where again the iteration is initiated by 
fa = 0.111n(100 GeV). 

In both cycles the iterated values are arguments of a logarithmic function; 
thus the iteration procedure is very fast. After the first iteration, an accuracy 
of about 0.1% has been achieved for energies in the range 80^-150 GeV, while 
a second iteration is needed to obtain the same precision for the other beam 
energies. In Fig. [7J the energy linearity, defined as the ratio between the mean 
reconstructed energy and the beam energy, is compared, after a first iteration, 
to the linearity obtained after iterating to a e = 0.1% accuracy, showing a good 
agreement. For this reason, the suggested algorithm of the energy reconstruc- 
tion can be used for the fast energy reconstruction in a first level trigger. 

Fig. [7] also demonstrates the correctness of the mean energy reconstruction. 
The mean value of E/E beam is equal to (99.5 ± 0.3)% and the spread is ±1%, 
except for the point at 10 GeV. However, as noted in [|TT|], result at 10 GeV is 
strongly dependent on the effective capability to remove events with interac- 
tions in the dead material upstream and to separate the real pion contribution 
from the muon contamination. 

Fig. H shows the pion energy spectra reconstructed with the e/h method 
proposed in this paper for different beam energies. The mean and a values 
of these distributions are extracted with Gaussian fits over ±2a range and are 
reported in Table [I] together with the fractional energy resolution. 

Fig. |9| shows the comparison of the linearity as a function of the beam energy 
for the e/h method and for the cells weighting method [pPp. Comparable quality 
of the linearity is observed for these two methods. 

Fig. [T^ shows the fractional energy resolutions (a/E) as a function of 1/ \fE 
obtained by three methods: the e/h method (black circles, also presented on the 
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Table [T]), the benchmark method [ JTT] ] (crosses), and the cells weighting method 
[ |TT| ] (open circles). The energy resolutions for the e/h method are comparable 
with the benchmark method and only 30% worse than for the cells weighting 
method. A fit to the data points gives the fractional energy resolution for the 
e/h method obtained using the iteration procedure with e = 0.1%, 

a/E = [(58 ± 3)%/VE + (2.5 ± 0.3)%] © (1.7 ± 0.2)/ E (9) 
for the e/h method using the first approximation, 

o/E = [(56 ± 3)%/VE + (2.7 ± 0.3)%] © (1.8 ± 0.2)/ E, (10) 
for the benchmark method, 

a/E = [(60 ± 3)%/VE + (1.8 ± 0.2)%] © (2.0 ± 0.1) /E, (11) 
and, for the cells weighting method, 

a/E = [(42 ± 2)%/y/E + (1.8 ± 0.1)%] © (1.8 ± 0.1)/ E, (12) 

where E is in GeV and the symbol © indicates a sum in quadrature. The sam- 
pling term is consistent between the e/h method and the benchmark method 
and is smaller by a factor of 1 .5 for the cells weighting method. The constant 
term is the same for the benchmark method and the cells weighting method and 
is larger by (0.7 ± 0.3)% for the e/h method. The noise term of about 1.8 GeV 
coincide for all four cases within errors that reflect its origin in electronic noise. 
Note, that from the pedestal trigger data the total noise for the two calorimeters 
was estimated to be about 1 .4 GeV. 



4 Hadronic Shower Development 

The e/h method for energy reconstruction has been used to study the energy 
depositions, Ei, in each longitudinal calorimeter sampling. Table lists (and 
Fig. [n] shows) the differential mean energy depositions (AE/Az)i = Ei/Azi 
as a function of the longitudinal coordinate z for energies from 10 to 300 GeV, 
with z expressed in interaction length units. 

A well known parameterization of the longitudinal hadronic shower devel- 



opment from the shower origin is suggested in Ref. [30]: 



N L f')- 1 e -»* + (1 _ w) (I.)- 1 e -«*7 > , (13) 



dz I \X J VAj 
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where N is the normalization factor, and a, b, d, to are parameters (a = 
0.6165 + 0.3183 -\nE,b = 0.2198, d = 0.9099 - 0.0237 -\nE,uj = 0.4634). 
In this parameterization, the origin of the z coordinate coincides with shower 
origin, while our data are from the calorimeter face and, due to insufficient lon- 
gitudinal segmentation, the shower origin can not be inferred to an adequate 
precision. Therefore, an analytical representation of the hadronic shower lon- 
gitudinal development from the calorimeter face has been used pT|]: 



1 F 1 (l,a+l,{d-l)^)\, (14) 



where \Fi(a, (3, z) is the confluent hypergeometric function. Note that the for- 
mula ( |14| ) is given for a calorimeter characterized by its X and A/. In the 
combined setup, the values of X , A/ and the e/h ratios are different for elec- 
tromagnetic and hadronic compartments. So, the use of formula ([14]) is not 
straightforward for the description of the hadronic shower longitudinal profiles. 



To overcome this problem, Ref. Q320 suggests an algorithm to combine 



the electromagnetic calorimeter (em) and hadronic calorimeter (had) curves 
of the differential longitudinal energy deposition dE/dz. At first, the mean 
hadronic shower develops according eq. ( |T4] ) in the electromagnetic calorimeter 
to the boundary value z em which corresponds to a certain integrated measured 
energy E em (z em ). Then, using the corresponding integrated hadronic curve, 
E(z) = f*(dE/dz)dz, the point Zhad is found from the equation Ehad(zhad) = 
E e m(z em ) + Ed m . From this point a shower continues to develop in the hadronic 
calorimeter. In principle, instead of the measured value of E em one can use the 
calculated value of E em = f* ern (dE / dz)dz obtained from the integrated elec- 
tromagnetic curve. The combined curves have been obtained in this manner. 



Fig. [TT| shows the differential energy depositions (AE/Az)i = Eij Azt as 
a function of the longitudinal coordinate z in units of A^ for the energy from 
10 to 300 GeV and a comparison with the combined curves for the longitudinal 
hadronic shower profiles (dashed lines). The level of agreement was estimated 
using the function x 2 where, following Ref. [^Up, the variances of the energy 
depositions are taken to be equal to the depositions themselves. A significant 
disagreement (P(x 2 ) < 0.1%) has been observed between the experimental 
data and the combined curves in the region of the LAr calorimeter, especially 
at low energies. 

We attempted to improve the description and to include such essential fea- 
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ture of a calorimeter as the e/h ratio. Several modifications and adjustments 
of some parameters of the parameterization ( |14| ) have been tried. The con- 
clusion is that replacing the two parameters b and u in the formula dl4| ) with 
b = 0.22 ■ {e/h) cal /(e/h)' cal andu; = 0.6 ■ {e/n) cal / (e/-R)' cal results in a reason- 
able description of the experimental data. Here the values of the (e/h)' cal ratios 
are (e/h)' em ~ 1.1 and (e/h)' had « 1.3 which correspond to the data used for the 
Bock et al. parameterization [|3~0|]. The (e/7r)' ca/ are calculated using formulas 
i)and(|). 

In Fig. [12] the experimental differential longitudinal energy depositions and 
the results of the description by the modified parameterization (solid lines) are 
compared. There is a reasonable agreement (the probability of description is 
more than 5%) between the experimental data and the curves. Note, that pre- 
vious comparisons between Monte-Carlo and data have shown that FLUKA 
describes well the longitudinal shape of hadronic showers JTT| ] • 

The obtained parameterization has some additional applications. For exam- 
ple, this formula may be used for an estimate of the energy deposition in various 
parts of a combined calorimeter. This is demonstrated in Fig. [T3]in which the 
measured and calculated relative values of the energy deposition in the LAr and 
Tile calorimeters are presented. The errors of the calculated values presented 
in this figure reflect the uncertainties of the parameterization ([14]). The relative 
energy deposition in the LAr calorimeter decreases from about 50% at 10 GeV 
to 30% at 300 GeV. Conversely, the fraction in the Tile calorimeter increases as 
the energy increases. 



5 Conclusions 

Hadron energy reconstruction for the ATLAS barrel prototype combined calo- 
rimeter has been carried out in the framework of the non-parametrical method. 
The non-parametrical method of the energy reconstruction for a combined calo- 
rimeter uses only the e/h ratios and the electron calibration constants, without 
requiring the determination of other parameters by a minimization technique. 
Thus, it can be used for the fast energy reconstruction in a first level trigger. 
The value of the e/h ratio obtained for the electromagnetic compartment of 
the combined calorimeter is 1.74 ± 0.04 and agrees with the prediction that 
e/h > 1.66 for this calorimeter. The ability to reconstruct the mean values 
of particle energies (for energies larger than 10 GeV) within ±1% has been 
demonstrated. The obtained fractional energy resolution is [(58 ± 3)%/yE + 
(2.5 ± 0.3)%] © (1.7 ± 0.2) /E. The results of the study of the longitudinal 
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hadronic shower development have also been presented. 
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Table 1 : Mean reconstructed energy, energy resolution and fractional energy 
resolution for the various beam energies. 



Ebeam 


E (GeV) 


a (GeV) 


a/E (%) 


10* GeV 


9.30 ±0.07 


2.53 ±0.05 


27.20 ±0.58 


20* GeV 


19.44 ±0.06 


3.41 ±0.06 


17.54 ±0.31 


40 GeV 


39.62 ±0.11 


5.06 ±0.08 


12.77 ±0.21 


50 GeV 


49.85 ±0.13 


5.69 ±0.13 


11.41 ±0.26 


80 GeV 


79.45 ±0.16 


7.14 ±0.14 


8.99 ±0.18 


100 GeV 


99.10 ±0.17 


8.40 ±0.16 


8.48 ±0.16 


150 GeV 


150.52 ±0.19 


11.20 ±0.18 


7.44 ±0.12 


300 GeV 


298.23 ±0.37 


17.59 ±0.33 


5.90 ±0.11 



*The measured value of the 
*The measured value of the 



beam energy is 
beam energy is 



9.81 GeV. 
19.8 GeV. 
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Table 2: The differential mean energy depositions AE/Az (GeV/X n ) as a 
function of the longitudinal coordinate z for the various beam energies. 



N 


z 


Ebeam (GeV) 


depth 


(A.) 


10 


20 


40 


50 


1 


0.294 


5.45 ±0.08 


8.58 ±0.16 


14.3 ±0.2 


16.6 ±0.4 


2 


0.681 


4.70 ±0.08 


9.10 ±0.15 


16.7 ±0.2 


20.8 ±0.3 


3 


1.026 


2.66 ±0.06 


5.55 ±0.11 


11.1 ±0.2 


13.6 ±0.2 


4 


2.06 


1.93 ±0.03 


4.35 ±0.06 


8.99 ±0.08 


11.0 ±0.1 


5 


3.47 


0.87 ±0.02 


2.13 ±0.04 


5.29 ±0.06 


6.15 ±0.10 


6 


5.28 


0.18 ±0.01 


0.57 ±0.02 


1.50 ±0.03 


2.07 ±0.05 


7 


7.50 


0.025 ±0.003 


0.11 ±0.01 


0.32 ±0.01 


0.49 ± 0.02 


N 


z 


£ 6eam (GeV) 


depth 


(K) 


80 


100 


150 


300 


1 


0.294 


22.6 ±0.6 


28.4 ±0.6 


36.3 ±0.7 


61.3 ± 1.5 


2 


0.681 


30.4 ±0.4 


37.6 ±0.5 


53.5 ±0.8 


97.9 ± 1.7 


3 


1.026 


20.3 ±0.3 


25.7 ±0.4 


37.2 ±0.6 


68.9 ± 1.2 


4 


2.06 


18.0 ±0.1 


22.4 ±0.2 


33.9 ±0.3 


64.8 ±0.7 


5 


3.47 


11.9 ±0.1 


14.6 ±0.2 


23.3 ±0.2 


49.0 ±0.5 


6 


5.28 


3.66 ±0.06 


4.57 ±0.08 


8.18 ±0.13 


18.6 ±0.3 


7 


7.50 


0.86 ±0.03 


1.10 ±0.04 


2.04 ±0.06 


5.54 ±0.15 
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Figure 1: Three-dimensional cutaway view of the ATLAS calorimeters. 
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Figure 2: Schematic layout of the experimental setup for the combined LAr 
and Tile calorimeters run (side view). The S3 and S4 are scintillation counters, 
the BC3 is a beam proportional chamber, and the midsampler and the "muon 
wall" are scintillation hodoscopes. 
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(e/7i) LAr (e/ji) LAr 



Figure 4: The distributions of the (e/n) LAr ratio for beam energies of 20 and 
50 GeV (top row, left to right), and beam energies of 100 and 300 GeV (bottom 
row, left to right). 
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Figure 5: The mean value of the (e/71") lat ratio as a function of the beam energy. 
The curve is the result of a fit of equations (Q) and (Q). 
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Figure 6: The {eji^LAr ratios as a function of the beam energy, for e/h method 
(black circles) and for weighting method (open circles for |Q] and open squares 
for tlT0|]). The lines are the result of a fit of equations @ and (f|) with free e/h 
parameter and k = 0.11: solid line is for our data, dashed line is for the [^] data 
and dash-doted line is for the [ [TO] ] data. 



26 



E 

CO 
<U 
Si 

UJ 
UJ 



1.02 



0.98 



0.96 



0.94 




100 



200 300 
E b eam (GeV) 



Figure 7: Energy linearity as a function of the beam energy for the e/h method 
obtained using the iteration procedure with e = 0.1% (black circles) and with 
the first approximation (open circles). 
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Figure 8: The energy distributions for beam energies of 20 and 50 GeV (top 
row, left to right), and beam energies of 100 and 300 GeV (bottom row, left to 
right). 
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Figure 9: Energy linearity as a function of the beam energy for the e/h method 
(black circles) and the cells weighting method (open circles). 
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Figure 10: The energy resolutions obtained with the e/h method (black cir- 
cles), the benchmark method (crosses) and the cells weighting method (circles). 
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Figure 1 1 : The comparison between the experimental differential mean energy 
depositions at 10 GeV (crosses), 20 GeV (black top triangles), 40 GeV (open 
squares), 50 GeV (black squares), 80 GeV (open circles), 100 GeV (black cir- 
cles), 150 GeV (stars), 300 GeV (black bottom triangles) and the calculated 
curves ( |T4"| ) as a function of the longitudinal coordinate z in units A^. The er- 
rors on the data points are the errors of the mean values and there are within 
symbols size. 
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Figure 12: The experimental differential mean longitudinal energy depositions 
at 10 GeV (crosses), 20 GeV (black top triangles), 40 GeV (open squares), 50 
GeV (black squares), 80 GeV (open circles), 100 GeV (black circles), 150 GeV 
(stars), 300 GeV (black bottom triangles) energies as a function of the longitu- 
dinal coordinate z in units for the combined calorimeter and the results of 
the description by the modified parameterization for the equation (J14T). 
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Figure 13: Energy deposition (percentage) in the LAr and Tile calorimeters at 
different beam energies. The circles (squares) are the measured energy depo- 
sitions in the LAr (Tile) calorimeter, the diamonds (crosses) are the calculated 
energy depositions. 
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